Dynamic clustering in active colloidal suspensions with chemical signaling 
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In this paper, we explore experimentally the phase behavior of a dense active suspension of self- 
propelled colloids. In addition to a solid-like and a gas-like phase observed for high and low densities, 
a novel cluster phase is reported at intermediate densities. This takes the form of a stationary 
assembly of dense aggregates, with an average size which grows with activity as a linear function 
of the self-propelling velocity. While different possible scenarii can be considered to account for 
these observations - such as a generic velocity weakening instability recently put forward -, we show 
that the experimental results are reproduced by a chemotactic aggregation mechanism, originally 
introduced to account for bacterial aggregation, and accounting here for diffusiophoretic chemical 
interaction between colloidal swimmers. 



PACS numbers: 

Active systems refer generically to collections of par- 
ticules which consume energy at the individual scale in 
order to provide self-propelled motion. In assembly these 
systems usually exhibit a wide variety of collective behav- 
iors, structures and patterns, which depart strongly from 
the classical equilibrium expectations [TMTQ] . In particu- 
lar, there are many observations in nature of cluster-like 
phases in very different systems: flocks, schools, swarms 
of fishes, insects or bacteria [TTHT6] . However this appar- 
ent analogy, occuring over a broad range of scales, hides 
many different mechanisms in the propulsion and inter- 
actions. In order to disentangle the universal from the 
specific behaviors of those complex phases, a systematic 
experimental exploration of artificial active particles sys- 
tems at high density is needed. Some experiments have 
been carried out on self-propelled walkers at high den- 
sities [17 but to our knowledge, active particles at the 
colloidal scale - involving natural brownian noise and 
solvent induced interactions - were only studied at low 
densities [IBM]. 

In this paper we explore experimentally the behavior 
of a two dimensional dense active suspension of artifi- 
cial self-propelled colloids (SPC). We first characterize 
the phase behavior of this active system under an exter- 
nal (gravity) field, from the dilute gas to a dense solid- 
like phase. A key observation is the emergence of dy- 
namic clustering at intermediate densities. Clusters of 
SPC form naturally in the system, with an average size 
which grows with activity, in direct proportionality to the 
propelling velocity of an individual SPC. Several scenarii 
are discussed in order to rationalize these experimental 
results, suggesting in particular a possible chemotactic 
aggregation mechanism. 

Experimental description - The active particles are 
home-made spherical gold colloids of radius a c± 1/im 
half covered with platinum [22] . In presence of hydrogen 
peroxide the particles self-propel consuming H2O2 under 
a self-phoretic motion (a combination of diffusiophoresis 
and self-electrophoresis [20j [23] ) . We note that one key 
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FIG. 1: Top: (left) Schematic representation of the exper- 
imental system; (right) normalized surface density profiles 
p/pmax as a function of the position z along the cell (with 
(blue) and without (red) H2O2). Solid lines are fits with a 
tangent hyperbolic function. Bottom: picture of the 2D sedi- 
mented particles, left: passive colloids in water; right: active 
colloids in a solution of 0.1% of H2O2. The effective gravity, 
g.sinO ~ 2.10 -2 m.s -2 and is indicated by the arrow. The 
pictures were taken with a 20 x objective. 



aspect here is that the system does self-propel - with a 
typical speed V ~ 3/im.s -1 - even with very low concen- 
trations of hydrogen peroxide. Typical volumetric con- 
centrations of H2O2 are in the range 0.01-0.1% - avoid- 
ing the formation of O2 bubbles even at high particles 
concentrations, which is a key issue for alternative active 
colloidal systems of same family-. The Peclet number 
for those active particles is typically Pe = V.cl/Dq ~ 13 
where Do is the bare diffusion coefficient of the colloids. 

We observe the system with an inverted optical micro- 
scope both in reflection and transmission modes. Videos 
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were taken with a Hammatsu Orca-ER camera at a fram- 
erate varying from 1 to 9 Hz. A Matlab routine al- 
lows particle detection and tracking. We confine those 
particles using a reduced gravity field: we let the col- 
loids sediment in a cell slightly tilted with an angle 
6 ~ 2.10 _3 rad. As the colloids are heavy (with mass 
density ~ 19 g.cm -3 ) they quickly sediment at the bot- 
tom of the cell and form a dense two dimensional layer 
(see Fig. [T]). As soon as we introduce hydrogen perox- 
ide the particles become active and self-propel, leading 
to strong change in the sedimentation profile. Profiles 
relax after a few minutes, allowing then to observe and 
characterize the system in a stable and stationary state 
over the whole duration of the experiment - typically a 
few hours-. 

Sedimentation profile and phase behavior - As observed 
in Fig. [TJ the sedimentation behavior of SPC is strongly 
affected by the amount of activity. Non-active SPC parti- 
cles suspended in water exhibit as expected a solid phase 
at the bottom, above which a very low density gas phase 
(at room temperature T am b) can be distinguished, see 
Fig. [l]-c (left). When SPC are made active by adding 
0.1% H2O2, Fig. [l]-c (right), a solid phase remains, but 
the top gas phase spreads to much higher heights, with 
an exponential decay of the density profile which can be 
accounted for by a high effective temperature T e ff ~ 50T, 
in agreement with the estimate based on the Peclet num- 
ber, see [20]. Furthermore, as can be observed in the ex- 
perimental snapshots in Fig.[l]-c, particles tend to cluster 
for intermediate densities, a behavior which we will ex- 
haustively explore below. 

The difference between passive and active particles can 
also be clearly seen on the surface density profile p of the 
colloids, Fig. [l]-b. These density profiles are obtained af- 
ter a binarization of the 500 movie frames (taken at lfps) 
and averaged over the horizontal direction. The experi- 
mental profiles exhibit an enlargement of the "interme- 
diate zone" between the solid and gas phases when the 
particles are made active. We finally note that the situ- 
ation is reversible which means that when we replace the 
hydrogen peroxide solution by water, this intermediate 
zone disappears to recover the non-active sedimentation 
profile. 

We then measure the structure factors of the system 
for various zones in the sedimentation profile, Fig. [2] The 
structure factor is defined as S(k) = -^{pkP-k) where N 
is the number of particles in the region of interest and 
pk is the Fourier transform of the instantaneous num- 
ber density of SPC. For the SPC without activity - in 
water -, we report S(k) in the bottom solid phase only. 
As shown in Fig. [2] (dotted line), it exhibits a strong or- 
dering characteristic of a solid phase. It is however not 
perfectly crystalline due to the size polydispersity of the 
colloids (typically around 10 %). For the active SPC, we 
have measured S(k) in three different zones arbitrarily 
defined as: (i) a solid phase - for p/pmax > 0.8 -, (ii) 
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FIG. 2: Structure factors S(k): Non-active SPC in water in 
the solid phase, dotted line; Active SPC (in H 2 2 = 0.1%) 
for: solid phase (green), intermediate zone (blue), gas phase 
(red). Inset: same curve in logarithmic scales, zooming in the 
small k behavior. 



an "intermediate zone" - for 0.05 < p/pmax < 0-8 ~ and 
(hi) a gas phase - for p/p 

max 

< 0.05 [37]; pmax is the 
maximum surface density (typically of order of the ran- 
dom close packing value). As seen in Fig. [2] while a high 
ordering remains in the solid phase of SPC, the ampli- 
tude of the first peak decreases strongly and activity is 
observed to destabilize the crystalline phase, in line with 
recent simulations [10] . For the low density phase, a sig- 
nature of ordering in the gas phase is also observed (up 
to 3 peaks can be identified). The intermediate "clus- 
ter" phase is particularly interesting: it exhibits a strong 
ordering, but also shows a strong increase of the struc- 
ture factors when k —> 0, Fig. [2j This is a signature of 
a strong compressibility of the system, reminding of a 
critical behavior and large fluctuations. 

Emergence of dynamic clustering - A key observation 
in the snapshots, Fig.[3]and movies [36] , is the appearance 
of clustering for active SPC in the regime of intermediate 
densities, corresponding typically to a surface fraction <p 
in the range 3% — 50%. In order to better characterize 
this 'cluster phase', we use a non-tilted cell and start with 
a homogeneous suspension of particles in water, with sur- 
face fraction ~5%. Upon addition of H2O2, we observe 
after a few seconds the appearance of clusters coexisting 
with a gas phase, see Fig. [3] Those clusters are not immo- 
bile but exhibit a slow, random-like, dynamics, with typ- 
ical speeds of 0.1/im.s -1 to be compared to V ~ 3pm. s -1 
the typical velocity of the SPC. Another striking point is 
that the clusters are dynamic, with particles going in and 
out over time, and clusters merging and dissociating, see 
Fig [3] and movie [36] . 

More quantitatively we measured the average size of 
cluster versus the activity of the system, i. e. the amount 
of fuel added, for volumetric concentration of hydrogen 
peroxide varying from 0.01% to 0.1%. At each H2O2 
concentration we measure the activity by tracking the 
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FIG. 3: Top: Dynamic clusters of SPC observed in a hor- 
izontal geometry. The surface fraction of SPC is 5% and 
H2O2 concentration is 0.1%. Images obtained in transmission 
mode with a 63 x objective. Top: snapshot of the clusters; see 
[36] for a movie. Bottom: chronophotography illustrating the 
dynamics inside one cluster at t = 0, 0.4 and 0.9 s. Marked 
particles show a dynamic exchange of SPC between the gas 
phase and the cluster, as well as internal reorganization of the 
cluster. 

individual SPC in the gas phase and extracting the mean 
velocity V. We count the number of colloids forming each 
cluster in the stationary state. A cluster is defined as an 
assembly of at least three colloids that keeps the same 
configuration for at least one second. For each hydrogen 
peroxide concentration, we count more than 150 different 
clusters. As shown in Fig.[4j we measure that the average 
cluster size N* is a linear increasing function of the mean 
velocity V of the individual SPC in the gas phase. This 
is a counter-intuitive result, as it shows that clustering 
is enhanced for larger effective temperature of SPC (as 
T eff ~F 2 [20]). 

Tentative scenarii for the dynamic clustering - The 
theoretical litterature offers various predictions for clus- 
tering in active matter, which involve a variety of ingre- 
dients and mechanisms [IH3j [5j [6] . We thus try to dis- 
criminate between various possible scenarii allowing to 
rationalize the emergence of clustering in our active SPC, 
especially for surface fraction as low as a few percents. It 
is accordingly important to first recall some facts about 
our SPC. Our SPC are janus spherical particles, which 
leads only to a negligeable anisotropy in the steric in- 
teraction between particles, and would a priori discard 
modelization of the observed behavior in terms of cou- 
plings between density and short-range orientational or- 
dering. A second important point is that self-propulsion 
takes its origin in our case in the self-phoretic motion of 
the colloids. This leads to a perturbation of the hydro- 
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FIG. 4: Mean cluster size vs the average velocity of the parti- 
cles in the gas phase, measured in the experiments (with SPC 
surface fraction 5%). The point at V = is measured for zero 
activity. The dashed line is a linear fit, N(V) = 1.6V + 1.4 
(with V in /im.s -1 ). 

dynamic velocity field around a SPC decaying as ~ 1/r 3 
[24] , in contrast to pusher-puller systems which behave as 
a force-dipole with a longer 1/r 2 decay. As for phoretic 
motion in general, hydrodynamic interactions are accord- 
ingly not expected to play a dominant role in our system. 
In contrast, chemical interactions are susceptible to act as 
a novel ingredient for our SPC, since, due to the chemical 
origin of the motion, the consumption of fuel and asso- 
ciated production of 'waste' species leads to long range 
chemical gradient around SPC [23 - with concentration 
fields decaying here like 1/r [25 -. Such gradients are 
expected to induce diffusiophoretic interactions between 
particles, as we discuss below. 

Among the exhaustive zoology of aggregation behavior 
discussed in theoretical litterature, many of the interpre- 
tations are difficult to reconcile with the present obser- 
vations: 

- Flocking behavior, i.e. group of orientationally ordered 
particles moving in the same direction [TJ[5]: in our case, 
SPC clusters do not exhibit directed motion. This dis- 
cards the interpretation in terms of a nematic ordering 
of the SPC, in line with the expected weak anisotropy of 
their steric interaction. 

- Adhesive clustering: one may a priori not discard in 
an obvious way a mechanism based on a simple adhesion 
between particles, by analogy to the structure of simple 
adhesive colloids [26 . However the observed SPC clus- 
ters are reversible and strongly kinetic by nature. This 
adhesive scenario is also difficult to reconcile with the ob- 
servation that the size of the SPC clusters increases with 
the activity of the particles, Fig. [4] 

- Surface generated slip flows: another mechanism could 
be the surface-generated attraction between particles in- 
duced by slip flows, as previously described for ther- 
mophoresis, Ref. [27] [28]. Such a scenario was shown in 
[27] to induce two dimensional crystalization of colloids 
at the surfaces. This is however quite in disagreement 
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with the present observation of kinetically alive clusters, 
reaching quickly a steady state, stable over hours. 

While the above scenarii do not provide a convincing 
framework to understand the present observations, we 
highlight two alternative views which account better for 
our experimental results. 

- First, in a recent numerical study, Fily and Marchetti 
[29] have shown that self-propelled particles with no 
alignment mechanism can exhibit an intrinsic cluster- 
ing instability, though at higher density. They interprete 
their results in terms of a generic instability taking its 
origin in the density-weakening dependence of the par- 
ticles velocities. Many features of our experiments are 
reproduced in their simulations and model, such as the 
stationary clustering state for intermediate densities and 
the corresponding strong increase of the structure factor 
for low fc, in direct line with our Fig. [2] However at this 
stage no prediction for the average cluster size N(V) is 
proposed which could account for the present linear de- 
pendence N(V) ~ V; 

- A second possibility can be also drawn, which is directly 
related to the chemical interactions mentionned above. 
The role of chemical sensing has already been studied 
in the context of patterns and clusters formation of bac- 
teria [13] HU [30j [31] . In a pioneering work Keller and 
Segel (KS) rationalized such behaviors on the basis of a 
mean-field description taking into account the diffusion 
of bacteria, a drift induced by chemical sensing and the 
production and diffusion of a chemo- attract ant [30] l3T]: 
its extension to thermotaxic SPC was proposed recently 
[32] . In our context a direct analogy may be made: the 
chemotactic effect is played by the diffusiophoretic mo- 
tion [33 induced by the gradient of species involved in 
the redox chemical reaction that propels the colloids [23] . 
We observed independently that the consummed H2O2 
act as a diffusiophoretic repellent to colloids, and its re- 
action induced depletion is expected to contribute to an 
overal chemo- attractive drift between SPC. Within the 
KS model, the spatio-temporal dynamics of the active 
particle population can be described at a mean field level 
in terms of a (here, 2D) particle density p and a global 
chemoattractant field c: 

d t p =D p V 2 p-V(fipVc); 

d t c =D c V 2 c + ap. (1) 

Here D p is the effective diffusion coefficient of the SPC, 
D c is the 'chemoattractant' diffusion coefficient, fi the 
diffusiophoretic mobility, a the chemical rate of the pow- 
ering chemical reaction occuring at the surface of each 
colloid. 

An interesting feature of the KS equations is that 
they exhibit singular solutions, leading to a 'chemotactic 
collapse' of the structure into a single or many dense 
aggregates (hence termed clumping) [31] • This phe- 
nomenon introduces a threshold 'Chandrasekhar' num- 
ber, N C1 above which the bacteria population clusterizes, 



while below it remains homogeneous. The expression for 
N c is (in 2D): A^ c = AD p D c / \iol. Furthermore, in the 
case of clumping - as discussed in [31] -, N* « N c is 
also expected to fix the typical size of a cluster. It also 
interesting to note that N* scales like the inverse of the 
self-trapping coupling constant g introduced by Tsori and 
de Gennes [34] , so that for a cluster g e # = g • N* ~ 1 is 
at the threshold for self-trapping and reduced mobility. 
Here, in order to connect this size to the experimentally 
measured velocity V of the SPC, we first note that V 
is itself a function of these parameters. Indeed the mo- 
tion is driven by the chemical reaction at the surface of 
the colloids and following [23] [25] , one expects typically 
V ~ fia/aDc, so that jia ~ V aD c . Furthermore, the 
effective diffusion coefficient of the particles scales like 
D p oc V 2 r r , with r r the rotational brownian time of the 
particle [20]. Altogether, this yields: 

Vt 

N* ~ — — ^ ^ Pe, (2) 

where Pe = V.cl/Dq is a Peclet number characterizing 
the SPC, which can be also interpreted in terms of a per- 
sistence number [35]. This prediction is in good agree- 
ment with the experimental result in Fig. [3j N(V) oc V, 
with furthermore a predicted prefactor given by an in- 
verse velocity r r /a ~ /im _1 .s, also in agreement with the 
experiments. 

This chemotactic scenario, relevant for chemically pow- 
ered SPC particles, proposes a consistent and predic- 
tive explanation for clustering. As a coarse-grained and 
mean-field description, KS cannot describe the detailed 
kinetics of the aggregate formation and distribution, 
which remain to be studied theoretically in more details. 
It would be also desirable to explore further the relative 
dynamic stability between a collection of cluster versus 
the relaxation towards a single cluster. Recent simula- 
tions of a chemotactic model have shown results similar 
to the one experimentally observed here with 'hot clus- 
ters' coexisting and fluctuating in time [35]. Altogether 
this framework provides interesting clues to be further 
explored theoretically in order to rationalize the cluster- 
ing mechanism observed in our work. 

To conclude, we have explored the high density phase 
behavior of suspension of active particles at the colloidal 
scale. We have observed the formation of a dynamical 
"cluster phase" with a typical cluster size increasing lin- 
early with a Peclet, or persistence, number characterizing 
the self-propelled colloids. Our results can be interpreted 
in the context a chemotactic aggregation scenario first 
introduced to explain clustering observed for bacterial 
population [13] [30]. It suggests that chemical interac- 
tions between SPC can mimic, on a purely physical basis, 
chemoattractivity and its consequences. 
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